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Jet tomography is the analysis of the attenuation pattern of high transverse momentum hadrons 
to determine certain line integral transforms of the density profile of the QCD matter produced in 
ultra-relativistic nuclear collisions. In this letter, we calculate the distortion of jet tomography due to 
multi-gluon fluctuations within the GLV radiative energy loss formalism. We find that fluctuations 
of the average gluon number, (N 9 ) ~ 3 for RHIC initial conditions, reduce the attenuation of pions 
by approximately a factor Z ~ 0.4 —0.5. Therefore the plasma density inferred from jet tomography 
without fluctuations must be enhanced by a factor 1/Z ~ 2. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 


I. INTRODUCTION 


The discovery of a factor of ~ 3 suppression of moderate pt < 4 GeV 7r°’s in central Au + Au reactions by 
PHENIX [ Q and large transverse asymmetries in non-central collisions for pt $ 5 GeV by STAR [ |2| have inspired 
several attempts [iilil to begin jet tomographic analysis of the matter density produced in ultra-relativistic 
nuclear reactions. Jet tomography is the QCD analog of conventional X-ray or positron tomography in that it exploits 
the attenuation of high energy jets produced in a variable density medium [ ||. The main source of attenuation of jets 
in QCD is induced gluon radiation due to multiple interactions in the medium. While the present data at moderate 
Pt are not conclusive, the tomographic analysis suggests that densities up to 100 times nuclear densities may have 
already been achieved. Soon very high statistics data out to pt ~ 10 — 20 GeV will be obtained, and it is important 
to refine the theory of jet tomography to take into account many sources of distortion of the attenuation pattern. In 
this letter we provide details of the calculation of the distortions due to gluon number fluctuations and compute the 
correction factors to the deduced densities. Such density renormalization factors Z ~ 0.4 — 0.5 have already been 
applied to calculations of particle spectra at RHIC [ Pi |l0|. 

Jet quenching probes the gluon density of the medium, characterized through the opacity parameter L/A g = 
I CTg(r)p s (r) dr. For the moderate opacities expected in nuclear collisions in the y/s ~ 200 AGeV range, it is 
convenient to calculate the induced radiation as a power series expansion in the opacity. In Ref. [ hi] we derived 
an analytic expression for the medium induced gluon radiation spectrum, p(x) = dN g /dx = p^Jx), in such a 

series form, where x is the light cone momentum fraction carried by the radiated gluon (see also [ [l^| and [0 for 
the relation of this approach to the asymptotic approach of Ref. [ |lj] ). 

The result for the n th order in opacity contribution to p(x) for a jet of energy E and color Casimir C R that is 
produced inside a medium with opacity L/A g can be expressed from [ [llj in the following form 
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The distribution of transverse momentum impulses is given by an effective dipole-like differential cross section that 
has an elastic component assumed to be given by a color screened Yukawa form [ |^] and a d-function component in 
the forward (jet) direction TqAn 2 ) 2 ~ We take p ~ gT ~ 0.5 GeV for RHIC initial conditions. 

One important advantage of the opacity formalism is that the effects of a finite temperature QCD plasmon frequency 
cutoff ujpi ~ gT /vl of the gluon radiation can be taken into account approximately via the kinematical cut-off both at 
small x < Xq = and at small k < p = iv m i n - We ignore the ~ a/ 3 difference between the plasma frequency and the 
screening mass since g ~ 2 and approximate both by ~ 0.5 GeV. The upper bound K 2 n . AX (x) = 4 E 2 min(x 2 , (1— x) 2 )—p 2 
in k 2 results from requiring the quenched jet as well as the radiated gluons have positive forward momenta. The 
fraction x c /x = p 2 L/(2xE) is a measure of the thickness of the medium to the gluon formation length. 

The radiation amplitudes are here denoted by A n ^ m = 2 (k— Q„) • (C„ im +i — C n , m ), where Q n = qi H-bq n , with 

Qo = 0, and C n , m = ^Vk log(k — Q„ + Q m _i) 2 . Destructive interference suppresses radiation with formation times 
greater than the thickness, L , of the medium. In addition to one power of x c /x in Eq.(l), higher order contributions 
are further suppressed by a formation factor from Eq. (116) of [ [hy 
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This simple analytic form arises for an exponential distribution, oc exp (zk — Zk- i)^rj> between adjacent scat¬ 
tering centers, Zk, in a plasma with mean thickness L at n th order in opacity. 

Gluon reabsorption from the medium reduces the radiation density p[x) at low x as shown in [ E- However, this 
effect is only important for jets with E less than a few GeV. Here we focus on higher energy jets. We consider only 
the static plasma geometry in Eq. (2) to simplify the numerical evaluation of p^ at higher orders. In Refs. [ 11-1 
we showed that the the main effect of 3+1D (Bjorken+transverse) expansion in the opacity expansion is to reduce 
the mean radiative energy loss, A E = ■ AE sta t. relative to the static approximation by a renormalization factor 

Z 3 +I = 2tq/L, where To is the formation time of the matter. The effective static opacity L/X = 5 that we use to fit 
the PHENIX data is relatively small because it value reflects the rapid dilution effects due to expansion. 

The numerically computed mean number of radiated gluons and the mean energy loss up to third order, 
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are shown in Figs. 1 and 2. The curves are calculated using p = 0.5 GeV, A g = 1 fm, and Cr = 3 and running a s . 
The leading first order gluon energy loss A E^\ is found to roughly follow the leading log expression A E^/E ~ 
^i^logjk f rom [ I®- The second order contribution reduces the first order result by a factor ~ 2 for E < 10 
GeV. By 40 GeV the second order correction is only 10%. For E > 40 both the second and third order correction 
are negligible, but also for E ~ 5 — 20 GeV the third order contribution largely cancels the second order one. The 
summed first through third orders shows that the induced fractional energy loss varies from 0.5 at E ~ 5 GeV down 
to ~ 0.3 at 20 GeV. Jet tomography at RHIC is dominated by quark fragmentation for pions with pt > 5 GeV with 
A E reduced by a factor Cf/Ca = 4/9. 



FIG. 1. A E/E is plotted versus E for opacity L/X g = 5. The three curves correspond 
to calculations up to 1 order (upper bound), 1+2 order (lower bound), and 1+2+3 
order (the actual final result). 

A similar pattern is seen in the energy dependence of the average number of gluons radiated in Fig. 2. The second 
and third order corrections largely cancel and the final average gluon number is only 2-3 for the kinematic range 
accessible at RHIC. Even with sensitivity ~ 2 on the kinematic cut-offs, ( N 9 ) is small. In fact as a function of energy 
the gluon number saturates for any given opacity L/\ g . The gluon radiative distributions are strongly peaked at 
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small x (which is consistent with the small x approximations used) and naturally need a lower cut-off generated by a 
characteristic jet energy independent mass scale to, i.e. x m i n ~ m/E. For any gluon radiative distribution that has a 
form p(x , E) « 1/a : c f(x/x c ) (with x c = typical energy scale/E) the mean number of gluons due to induced radiation 
(N 9 ) is approximately jet energy independent. As the energy of the jet and the radiative energy loss increase ( N 9 ) 
remains small but the radiated gluons become harder. 

We note that both A E and ( N 9 ) are much smaller than the recent estimates in [ 0 because our effective transport 
coefficient p 2 /\ g ~ 0.25 GeV 2 /fm, which is constrained by our fit below (see Fig. 4) to the PHENIX pion attenuation 
data [ |J, is 4 times smaller than the one considered for illustration in [ [bj|]. As noted before, expansion greatly 
reduces estimates based on the transport properties of the high density initial conditions. In addition, the finite 
plasmon frequency cut-off, ~ p, is the medium regulator of soft gluon number divergences. We emphasize that it is 
the finite medium dependent cut-off of the soft spectrum together with the relatively small transport coefficient and 
opacity of the plasma produced at RHIC energies that allows the opacity expansion to converge so rapidly even under 
the extreme conditions produced at RHIC. 



FIG. 2. Average number of gluons, ( N 9 (E )), with u> > p is plotted versus E for 
opacity L/\ g = 5. A set of curves corresponding to the curves in Fig. 1 is presented. 

II. FLUCTUATION SPECTRUM OF RADIATIVE ENERGY LOSS 

In the approximation that the fluctuations of the gluon number are uncorrelated, the spectrum of the total radiative 
energy loss fraction, e = u>i/E, can be expressed via a Poisson expansion P(e, E) = P) with Pi(e, E) = 

e~( N9 ) p{e, E) and 
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The form of this spectrum guarantees that the mean value is as in Fig. 1: 
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The above distribution differs considerably from that computed in [|l4j because ( N 9 (E )) is finite and small in our case 
and because p{x) = 0 for x < p/(2E) due to the plasma frequency cutoff that we impose on the low frequency modes 
in the medium. Therefore, we have explicitly a finite n = 0 (no radiation) contribution Po(e,E) = e~^ N 6(e). 
Assuming negligible kinematic correlations the numerical iteration of the recursion relation in Eq. (4) becomes very 
fast and is computed to high order (n < 25). However, by not enforcing that Xi < 1, there is a “leakage” error 
into the unphysical e > 1 range. We calculate this “leakage” error J) de P(e,E) and correct the normalization of 
P(e, E) in the the physical range e £ [0,1]. 

The resulting spectrum (without the delta function contribution) is shown in Fig. 3 for three different jet energies. 
The finite intercept at e = 1 provides a measure of the “leakage” error, which is acceptable in this case. The low 
frequency plasmon cut-off at ~ p, is clearly visible. In addition multi-gluon iterations in the probability distribution 
P(e,E) exhibit a slight oscillatory pattern in multiples of Xq. The high frequency random oscillations provide an 
indication of the accuracy of our Monte Carlo numerical integrations methods. 

The results indicate that P(e,E) is approximately constant from x 0 = p/(2E) up to a scale ~ x c = p 2 L/(2E). 
For x x c , P(e,E) decreases rather quickly. The (normalized to unity) probability distribution per gluon, p(x,E = 
40 GeV)/ ( N 9 ) is also shown for comparison. Multi-gluon fluctuations flatten the rapid small x rise of p even though 
Eqs. (||,||) dictate that the first moment of both p(x, E) and P(e, E) are the same. 



FIG. 3. Probability density of total fractional energy loss e = y~) uij/E for a gluon get 
with E=10,20, and 40 GeV traversing matter with opacity L/X g = 5. The numerical 
curves include orders 1+2+3 order in the opacity expansion for p(x). The n = 0 no 
radiation delta function contribution, Tb( e ), is not shown above. The low frequency 
plasma cut-off is at e = p/(2E). The probability density per gluon is also included. 


III. THE QUENCHING PATTERN OF tt° 

We apply the energy loss spectrum to calculate the quenched spectrum of hadrons by modifying the mean energy 
loss pQCD formulas from Refs. [ I! 110 0. We concentrate on mid-rapidity hadron production (y C m = 0). 
A jet of flavor c and transverse momentum p c produced in a hard PQCD scattering a + b — > c + dis attenuated 
prior to hadronization by the radiative energy loss to p* = p c ( 1 — e). This shifts the hadronic fragmentation fraction 
z c = Ph/Pc to z* = z c /( 1 - e). 

The invariant distribution of 7r° reduced by energy loss in central A + A collision is then given by 
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where Taa( 0) is the Glauber profile density in central collisions. The pion fragmentation function D^o/ c {z,Q 2 ) is 
taken from BKK [ |l8) . We take the GRV94 LO [ [H| structure functions for f a / P (x, Q 2 ) and include isospin dependence 
(Z protons and A — Z neutrons). Nuclear shadowing, intrinsic h? broadening and Cronin effect can be taken into 
account as in [ |2(j. |22|. ^3| . The interplay between the soft and hard components of hadron production studied in [ 

0.0 lead to modifications of the spectral shapes in the low px region and are neglected in this analysis. The factor 
z*/z c appears because of the in-medium modification of the fragmentation function [ Thus, the invariant cross 
section Eq. (j|) depends on the average opacity L/X g through the effect of P(e,p c ). 

We consider three different approximations to P(e,E): 

1. Use only the mean energy loss with P(e, E) « S(e — AE(E)/E) as in [ |], ||, ||, [|, |], [l7j 

2. Use the full fluctuating spectrum, P(e,E), from Eq. (|j) 

3. Use a renormalized average energy loss with P(e, E, Z) ss S(e — Z ■ A E(E)/E) 



FIG. 4. Quenching pattern for n° versus transverse momentum order by order in 
the opacity expansion. Curves labeled < APi+... > are calculated using the average 
energy loss up to the order indicated. The average opacity is taken to be L/\ g = 5. 

The data are from PHENIX [ 0 - Curves labeled Pi+...(e) average over the fluctuat¬ 
ing distributions illustrated in Fig. 3. The dot-dashed curves correspond to using a 
reduced (Z renormalized) first order average energy loss. 

The ratio, Raa(pt), compares the quenched to the unquenched 7r° distributions. In the case (1), the convergence of 
the opacity series using the mean energy shift to first and up to third order appears to be reasonably fast even though 
the second order correction is still uncomfortably large below px $ 10 GeV. Improved numerical methods need to be 
developed to enable summing higher order terms to verify our expectation that the summed results to third order are 
not significantly changed by higher order due to the additional l/(n + 1) and P„ jm factors in Eq. (1). 

We see from Fig. 4 that with even the modest value of the opacity L/X g = 5, the mean energy loss approximation over 
predicts the observed quenching by about a factor of two. Including the fluctuations in the Poisson approximation 
via P(e,E) leads to less energy loss by approximately a factor of two and brings the attenuation in line with the 
observed results. This renormalization of the effective energy loss can be inferred from the dot-dashed curves using 
approximation (3) above with Z ss 0.4 — 0.5. We conclude that the distortion of jet tomography due to gluon 
number fluctuations in the Poisson approximation can be well approximated by renormalizing the mean energy loss 
calculations by a factor Z ~ 0.5. 

While the px range of the available data is still too low to draw definitive conclusions, the effective static opacity 
with gluon fluctuation renormalization above corresponds from the results of [ [f| to an estimated initial gluon 
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rapidity density dN 9 /dy ~ 800 ± 100 and implies that the initial gluon density produced at RHIC may have reached 
p g (dN 9 /dy )/(tq ttR 1 2 ) ~ 20/fm 3 ~ lOOp^- 
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